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DETERMINATION OF TRAJECTORY AND ANGLES OF ATTACK 


OF A SCOUT HEAT-TRANSFER SPACECRAFT DURING 

REENTRY FLIGHT IN THE ATMOSPHERE 

By Robert J. Mayhue 
Langley Research Center 


SUMMARY 

A five- stage Scout reentry spacecraft experienced large-amplitude motions 
during the experimental heating data period as the result of disturbances from 
a fourth- stage motor failure. Downrange tracking radar failed to acquire the 
thrusting spacecraft at separation and tracking data for the experiment were 
limited to partial optical measurements from Bermuda. In order to evaluate the 
heat-transfer data obtained during the experiment, an analysis of the space- 
craft unknown reentry trajectory and angles of attack was made primarily on the 
basis of telemetered measurements of the spacecraft linear accelerations and 
angular velocities. Computational procedures established for the analysis 
individually employed both the flight measurements and the aerodynamic charac- 
teristics of the spacecraft from preflight wind-tunnel tests. 

The re suit s of the analysis indicated that the spacecraft entered the 
heating data period with a velocity of about 22 230 feet per second at an alti- 
tude of approximately 308 000 feet. These flight conditions were 3 800 feet 
per second lower in velocity, and at an altitude which was 51 000 feet higher 
than desired as the result of fourth-stage motor failure. The spacecraft maxi- 
mum angle of attack at the beginning of the heating data period was approxi- 
mately 40°, but damped to about 23 0 at the time the nose cap burned off. After 
the nose cap burned off, the resultant angle of attack was reduced further to 
a maximum of about 6° at maximum dynamic pressure, and indicated that the 
spacecraft was statically and dynamically stable throughout the experimental 
heating data period. 


INTRODUCTION 


As part of a general research program on reentry heating, the 
Research Center has conducted a free-flight test of a blunt-nose conical space- 
craft using the five-stage Scout booster vehicle. A malfunction of the fourth- 
stage rocket motor during the flight test caused premature separation of the 
spacecraft with large -amplitude coning motions. Tracking radar at Bermuda failed 




to acquire the spacecraft at separation, and tracking data for the experiment 
was limited to 4- seconds of optical measurements of the spacecraft angular 
position relative to the Bermuda optical tracking station. 

The purpose of this report is to describe the results obtained, and the 
methods employed, in the determination of the unknown trajectory and angles of 
attack of the spacecraft that were required for evaluation of the heat -transfer 
measurements. These results were based primarily on analysis of rate gyro and 
accelerometer measurements recorded during the experiment and were used for 
analysis of the heat-transfer measurements as reported in reference 1. 

The analysis of this report involved many unknown initial conditions and 
required the use of trial-and-error procedures in order to obtain results from 
the flight measurements which were consistent with the wind-tunnel aerodynamic 
characteristics of the spacecraft. The procedures and results that are pre- 
sented and discussed are as follows: 

(1) Integration of the rate gyro and accelerometer measurements, 

(2) Simulation of the experimental measurements using wind-tunnel aero- 
dynamics, and 

(5) Comparison of the forces and moments computed from the flight measure- 
ments with those obtained from wind-tunnel tests. 
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SYMBOLS 

longitudinal acceleration, g units 

longitudinal acceleration of center of gravity (high-range 
accelerometer), g units 

longitudinal acceleration of center of gravity (low-range 
accelerometer) , g units 

normal acceleration of center of gravity, g units 

resultant acceleration of center of gravity, g units 

transverse acceleration of center of gravity, g units 

total-force coefficients along X, Y, and Z body axes, 
respectively 

resultant-force coefficient in body-axis Y,Z plane 
resultant moment coefficient in body-axis Y,Z plane 
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force coefficients along Y and Z body axes due to pitching and 
yawing velocities, respectively 

moment coefficients along Y and Z body axes due to pitching and 
yawing velocities, respectively 

maximum cross-sectional reference diameter, 1.673 ft 

total-force components along X, Y, and Z body axes, respectively, 
lb 

aerodynamic-force components along X, Y, and Z body axes, 
respectively, lb 

gravity- force components along X, Y, and Z body axes, 
respectively, lb 

acceleration due to gravity at altitude, ft/sec2 
acceleration due to gravity at surface of earth, ft/sec2 
altitude, ft 

mass moment of inertias about X, Y, and Z body axes, respectively 

North latitude measured positive from equatorial plane, deg 
or rad 

Euler angle direction cosines between spacecraft body-axis 
system and local gravity-axis system 

Mach number 

aerodynamic moments about X, Y, and Z body axes, respectively, 
ft-lb 

mass, slugs 

free-stream dynamic pressure, lb/sq ft 

angular velocities about X, Y, and Z body axes, respectively, 
rad/sec or deg/sec 

spherical earth radius, 20.890150 X 10^ ft 
maximum cross-sectional reference area, 2.197 sq ft 
thrust , lb 
flight time, sec 

inertial velocity components along X, Y, and Z body axes, 
respectively 

earth- relative velocity, ft/sec 
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W weighty lb 

X g ,Yg,Zg orthogonal gravity-axis coordinate system with origin at center 

of gravity 

X,Y,Z orthogonal body-axis coordinate system with origin at center of 

gravity 

x C g longitudinal position of center of gravity measured from nose, 

ft 

x C p longitudinal position of center of pressure measured from nose, 

ft 

a angle of attack, deg 

3 angle of yaw, deg 

T| resultant angle of attack, deg 

7 earth- relative flight-path angle, rad or deg 

e earth-relative heading angle measured positive from East to 

South, rad or deg 

A West longitude measured negative from Greenwich meridian, 

deg or rad 

p free -stream density, slug/cu ft 

t aerodynamic phase angle, deg 

cu e angular velocity of rotating earth, 7.292115 x 10-5 rad/sec 

9,^,$ spacecraft attitude angles relative to gravity-axis system, 

rad or deg 

fi resultant spacecraft attitude angle in 9,t plane, rad or deg 

£ amplitude of coning motion, rad or deg 

Subscripts: 

o initial conditions 

max maximum condition 

min minimum condition 

b burnout 


A dot over a symbol denotes differentiation with respect to time. Bars 
indicate earth- relative values. 
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SPACECRAFT AND INSTRUMENTATION 


Spacecraft Description 

The spacecraft configuration consisted of a symmetrical blunt-nose, low- 
fineness-ratio, conical body with a 9° half-angle as illustrated in figure 1. 

A NOTS 100B 17-inch spherical rocket motor was incorporated in the spacecraft 
as the final thrusting stage of the Scout ST-8 booster system. An inconel nose 
cap calorimeter was installed to measure heat-transfer data at the beginning 
of reentry into the atmosphere. This cap subsequently burned off to expose a 
new configuration with a Teflon nose cap. Thermal protection for the side 
walls consisted of a Teflon skin bonded on a magnesium shell. 

The mass characteristics of the dynamically balanced spacecraft were 
obtained from preflight measurements and are listed for each fifth-stage flight 
event in the following table: 


Table I 


Flight event 


Weight 

lb 


Fifth- stage ignition 
Fifth- stage burnout 
Inconel nose cap off 


296.25 

155-45 

11 * 6.97 


> x cg> 


from base 
in. 


14.89 

20.33 

19.74 


> Iy = I Z , 

slug-ft2 slug-ft^ 


6.83 


3-93 

3.67 


Spacecraft Instrumentation 


Telemetered instrumentation was installed in the spacecraft to provide 
continuous measurements of linear acceleration and angular velocity along each 
body axis. Accelerometers were alined and fixed near the spacecraft center of 
gravity, and rate gyros were alined and fixed on the axes about which the mass 
moments of inertia were taken. The body-axis coordinate system and the sign 
convention used for the measurements along each axis are shown in figure 2. 
Ranges selected for the instrumentation are given in the following table: 


Table II 


Measurement 

Range 

Longitudinal acceleration: 

(a) High range, a L ^ D 

(b) Low range, a L .p 

Transverse acceleration, a«p . . 
Normal acceleration, ajj .... 

Rate of pitch, q 

Rate of yaw, r 

Rate of roll, p 

+lg to -6og 
0 to lOg 
+5g to -5g 
+5g to -5g 

+100 deg/ sec to -100 deg/ sec 
+100 deg/sec to -100 deg/sec 
1,200 deg/sec 












The displacements of the accelerometers from the spacecraft center of 
gravity during fifth-stage flight events were measured prior to launch and are 
given in the following table : 


Table III 


Measurement 

Displacement from 
center of gravity,* ft, of - 

AX 

AY 

AZ 


Fifth- stage 

ignition 


a L,T 


0 

-0.0427 

a L,D 


0 

.0427 

3-rp 

•5215 

-.1538 

0 

a N 

.5221 

0 1 

-.0056 

Fifth- stage burnout 

a L,T 

0.2067 

0 

-0.0427 

a L,D 

.2067 

0 

.0427 

&1JI 

.0682 

-•1538 

0 

a N 

.0687 

0 

-.0036 

Nose cap off 

a L,T 

0.2559 

0 

-0.0 427 

a L,D 

• 2559 

0 

.0427 

a T 

•1173 

-.1538 

0 

a N 

.1179 

0 

-.0036 


*Sign convention shown in figure 2. 


Ground Facilities 

The primary data acquisition was obtained from facilities located at 
Coopers Island in the Bermuda tracking network. Tracking equipment consisted 
of an FPS-l6 radar and ballistic streak and sequence cameras for optical cover- 
age. Radar acquisition of the spacecraft at separation from the fourth stage 
was to have been accomplished by means of skin track. 

Meteorological data were obtained from Coopers Island with sounding rocket 
launches conducted before and after the flight test. 
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Accuracy- 


Errors in the rate-gyro and accelerometer measurements were estimated to 
be ±2 percent of the instrumentation full-scale range and are listed in the 
following table: 


Table IV 


Measurement 

Error 

Longitudinal acceleration: 

(a) High range, a^p 

(b) Low range, a L T 

Transverse acceleration, ajp .'. 
Normal acceleration, ajj .... 

Rate of pitch, q 

Rate of yaw, r 

Rate of roll, p 

±1.22g 
±0 . 20g 
±0. lOg 
±0.10g 

±2.0 deg/sec 
±2.0 deg/sec 
±2^.00 deg/sec 


No attempt was made to estimate the accuracy of the radar tracking meas- 
urements, since radar failed to acquire the spacecraft at separation and con- 
tinued to track the burned-out fourth stage. Errors in the optical tracking 
angular measurements obtained during the experiment were believed to be negli- 
gible and were therefore used directly as a primary checkpoint for trajectory 
computations. 

Measurements of the free- stream density and speed of sound obtained from 
meteorological instruments up to an altitude of about 200 000 feet were esti- 
mated to be accurate within approximately ±5 percent. These quantities were 
extrapolated to the higher reentry altitudes where the effect of these errors 
became less significant. 


FLIGHT TEST 


The nominal trajectory and flight events predicted for the experiment are 
described in figure 3- The mission requirement was to obtain a reentry rela- 
tive velocity of about 2 6 000 feet per second at the end of fifth- stage 
thrusting. The heating data period was to begin at this time at an altitude of 
about 257 000 feet with a relative flight-path angle of approximately -14.0°. 

The actual trajectory and events that occurred during the experiment were 
as expected up through spinup and ignition of the fourth-stage rocket motor. 
About halfway through the fourth-stage thrusting period, however, a loss of 
thrust occurred and was accompanied by large -amplitude oscillations of the 
vehicle. The analysis of reference 2 showed that this event was caused by a 
blowout of the fourth-stage motor case which produced large moment disturbances 
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and loss of rocket-motor chamber pressure. The moments during thrust termina- 
tion caused the vehicle to diverge to a maximum angle of attack of about 58°, 
and loss of chamber pressure caused premature ignition of the spacecraft spher- 
ical rocket motor through a pressure-decay arming switch. 

Bermuda radar which had been beacon tracking the vehicle failed to acquire 
the thrusting spacecraft as it separated and accelerated away from the fourth 
stage. Tracking data obtained during the experiment were therefore limited to 
optical measurements from Bermuda which were obtained about 1 6 seconds after 
burnout of the spacecraft spherical rocket motor. These measurements did not 
provide specific trajectory coordinates, but were restricted to angular meas- 
urements of the spacecraft azimuth and elevation relative to the tracking sta- 
tion at 1- second intervals from 43 6 seconds to 44o seconds of flight time. 

Telemetered accelerometer and rate-gyro measurements were recorded during 
the reentry flight of the spacecraft and are presented in figures 4 to 8. Fig- 
ure 4 presents time histories of the spacecraft angular velocities from fifth- 
stage ignition through about 73 seconds of coasting flight. Ignition of the 
fifth-stage motor occurred at 377*85 seconds with burnout at about 420.85 sec- 
onds. Predicted times for these events were 393*8 seconds and 436.8 seconds, 
respectively. Cross plots of the rate-gyro measurements during fifth-stage 
thrusting and the beginning of coasting flight are presented in figure 5* Fig- 
ure 5(a) illustrates the large-amplitude circular coning motions, with jet 
damping effects that were experienced by the spacecraft during thrusting of the 
fifth-stage motor. After motor burnout, figure 5(b) shows that the spacecraft 
entered the heating data period with residual coning motions, and continued to 
oscillate with increasing angular velocities as indicated by figures 5(c) and 
5(d) . A sudden change in the amplitude of the angular velocities was observed 
at about 439 seconds in figure 5(d) which indicated burnoff of the inconel nose 
cap. This event was verified by the start of heating measurements on the 
exposed Teflon nose as reported in reference 1. Continuation of the cross 
plots of the rate gyro measurements through maximum dynamic pressure is shown 
in figures 5(c) to 5(m) • 

Figure 6 presents time histories of the spacecraft normal and transverse 
accelerations along the Y and Z body axes, as well as the resultant of these 
accelerations in the Y,Z plane. Cross plots of these measurements during 
coasting flight in the atmosphere through maximum acceleration are presented in 
figure 7* Longitudinal accelerations along the X body axis are presented in 
figure 8 from fifth- stage ignition through maximum acceleration of about 
-48. 6g units. All accelerometer data presented in figures 6 to 8 were cor- 
rected for errors due to displacement of the instrumentation from the space- 
craft center of gravity. 

It should be noted in figures 4 and 6 that the ranges of the rate gyros 
and accelerometers were exceeded and that the peak values were faired during 
the time intervals indicated. Fairings through data noise periods were also 
required at other time intervals as shown. 

Figures 9 and 10 present free-stream atmospheric data obtained from 
Bermuda radiosonde and sounding rocket measurements. A comparison of the 
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measured values with 1959 ARDC 1962 standard atmosphere is included. This com- 
parison indicated that the measured density became increasingly lower than 
standard through the altitude range that measurements were obtained. Extrapo- 
lation of the free-stream density to the higher altitudes, therefore, was based 
on an increasing percentage deviation as shown in figure 9 . The values for 
speed of sound presented in figure 10 , based on the free-stream temperature 
shown, were used to determine Mach number variations during all reentry trajec- 
tory computations. 

The time history of the spacecraft weight during thrusting of the NOTS 
100-B spherical rocket motor was computed from the longitudinal accelerometer 
measurements and predicted specific impulse of the motor. No drag or thrust 
misalinement being assumed, the computed propellant weight time history is pre- 
sented and compared with predicted nominal values in figure 11. Since compari- 
son with the predicted nominal indicated good agreement, these data were used 
to compute time histories of the spacecraft weight, center of gravity, and 
moments of inertia as presented in figure 12. Initial values for these quan- 
tities were obtained from preflight measurements of the spacecraft mass 
characteristics. 

Figure 13 presents the thrust time history of the NOTS 100-B motor as com- 
puted from longitudinal accelerometer measurements and the spacecraft weight 
variation in figure 12. Comparison of the computed thrust with predicted nom- 
inal values indicated that the thrust level during flight was slightly higher, 
and thrust tail-off occurred about 3 seconds earlier than predicted. 


METHOD OF ANALYSIS 


Comput at ional Procedure s 

A knowledge of the reentry trajectory and angles of attack was required 
for analysis of the primary heat-transfer measurements for about the first 
19 seconds of coasting flight. Flight times for this data period were from 
burnout of the fifth-stage motor at about 421 seconds to loss of the calorim- 
eter nose cap at approximately ^39 seconds. Trajectory and angle-of-attack 
computations to 460 seconds were also required for evaluation of ablation meas- 
urements through maximum dynamic pressure. 

Computation of the unknown trajectory and angles of attack of the space- 
craft by means of the rate gyro and accelerometer measurements presented a 
problem involving many unknown initial conditions. Extensive trial-and-error 
analysis using the two independent computational procedures was therefore 
required. The first computation, described in detail in reference 2 and 
appendix A of this report, was based on continuous integration of the rate-gyro 
and accelerometer measurements which provided increments of body-axis inertial 
velocities during the flight. The spacecraft aerodynamics were not required in 
this computation, and the aerodynamic forces and moments derived from these 
results were dependent on initial conditions. Constraints were imposed, 
therefore, which required that the aerodynamic characteristics derived from 
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integration of the rate-gyro and accelerometer measurements must he consistent 
with those from preflight wind-tunnel tests of reference 3- 

In order to monitor results from integration of the flight measurements, 
and to conduct initial-condition studies, a second independent computational 
procedure using wind-tunnel aerodynamics was established for the analysis. 

This computation was a six-degree-of-f reedom analog simulation of the flight 
measurements during coasting flight in the atmosphere and is described in 
appendix B. 


Initial Conditions 

If the Bermuda radar measurements of the spacecraft position coordinates 
at fifth- stage ignition were used, the equations of appendix A show that com- 
putation of the trajectory and angles of attack from rate-gyro and accelerom- 
eter measurements required six unknown initial conditions involving the body- 
axis inertial velocities u Q , v 0 , and w G and the Euler attitude angles 0 O , 

and 0 O • 


Since the pitch and yaw rate-gyro measurements near fifth- stage ignition 
were essentially sinusoidal, circular coning motions being indicated, the 
initial Euler angles were varied in the analysis by trial-and-error displace- 
ment of the coning motion from the velocity vector. The objective of this 
investigation was to determine the spacecraft thrusting attitude variations 
which would result in flight through the azimuth and elevation angles from 
Bermuda that were measured l6 to 20 seconds after burnout. In addition, tra- 
jectory increments as the result of thrusting off the velocity vector were 
required that would generate a dynamic -pres sure time history compatible with 
wind-tunnel longitudinal-force data as illustrated by the following equation: 



Determination of a trajectory which would result in agreement with Bermuda 
optical tracking measurements depended on the time histories of the spacecraft 
latitude, longitude, and altitude which varied according to the following 
equations : 


L = 


r T7 


"\ 


- L ° + J - 


V cos 7 sin e 
R 0 + h 


dt 


A = / V cos 7 cos . e dt 

° J (R 0 + h)cos L 


( 2 ) 


h = h Q + / V sin 7 dt 
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Equation (l) shows that for a given variation of velocity, the 
longitudinal- force coefficient derived from the computed dynamic pressure was 
dependent only on the density-altitude relationship specified by the flight- 
path-angle variation. The spacecraft position coordinates computed from equa- 
tions (2) , however, were dependent on the combined variations of heading and 
flight-path angles. Therefore, various coning-axis locations in the 0,t plane 
were assumed (see following sketch) until the interchange of heading and flight- 
path-angle increments during thrusting, as well as the velocity and altitude 
increments, yielded a reentry coasting trajectory in the atmosphere that 
agreed with optical and wind-tunnel measurements. 



The amplitude of the initial displaced coning motion was estimated by 
assuming that no external moments were present near ignition, and that the 
motions were centered about a total angular momentum vector fixed in space at 
the coning axis. The initial coning amplitude was then defined by the fol- 
lowing expression for the angle between the spacecraft spin axis and the total 
angular momentum vector: 


= Iy\Ao 2 + r o 2 
Po I X 


(3) 


The starting point on the coning circle at 0 O , and the tangent 
slope Q q at this point were chosen graphically in order to generate motions 
about the selected coning-axis orientation. The initial roll angle 0 O 
required to produce the initial slope fi 0 was computed from the following 
expressions: 


0 O = % cos 9*0 ' r o sin 0o 
r Q cos 0 O + % sin 0 O 


i|r = 


cos 0 
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where solving for 0 O gives: 


and 


0o 


tan"-*- 


r o • 

COS Qp ° 




cos 0 O 



(5) 


(6) 


The initial pitch and yaw rates were obtained from rate-gyro measurements, 
and the quadrant of 0 O was determined by the sign of the individual slopes 
0 O and \jr D in equations (4) . 

For a given set of initial Euler angles, the direction cosines for trans- 
formation of quantities in the gravity-axis system to the body-axis system were 
computed from equations (A9) in appendix A. The initial body-axis inertial 
velocities were then determined from the following equations: 


u o _ *g 0 l l + ^go* 2 " 

v o = Xg 0 m l + Y go m 2 - h 0 m 3 > 

w o = X go n 1 + Y gQ n 2 - h Q n 3 


(7) 


where the inertial velocities in the gravity-axis system were found from the 
equations : 

Y~ 

Xg = — + (Rq + h) Q u> e cos L 0 

u tan e 0 

Y = 4q s in ip 
tan 7 q 

= V Q sin y Q 

The spacecraft position coordinates X gQ , Y go , L 0 , A 0 and velocity- 
vector magnitude and orientation V D , i Q , 7 Q required for solution and 

integration of equations (8) were obtained from Bermuda radar measurements at 
fifth- stage ignition. 


> ( 8 ) 


12 


Preliminary integration results .- An example of the results obtained from 
preliminary trial-and-error integration of the flight measurements during 
thrusting is shown in figure l4. The initial conditions at fifth- stage igni- 
tion used for this computation resulted in divergence in the motions after 
burnout and were unacceptable, since corresponding analog simulation indicated 
that the spacecraft motions near burnout always precessed toward and eventually 
damped around the velocity vector. Additional integrations were made by graph- 
ically rotating and translating the previous 0,\|r motion patterns until a pos- 
sible solution of the spacecraft trajectory and motions during thrusting was 
obtained as shown in figure 15- Motions during thrusting are shown only near 
ignition and near burnout for clarity. The significant results of this inte- 
gration was the discrepancy between the precessional characteristics of the 
motions after burnout that were computed from the flight measurements and wind- 
tunnel aerodynamics as shown in figures l6(a) and l6(b) . Analog simulation 
results using wind-tunnel aerodynamics in figure l6(b) indicated that the oscil- 
lation amplitudes were larger than those obtained from integration of the 
flight measurements in figure l6(a) . From the results of figure 17; which pre- 
sents a comparison between the angles of attack computed from faired and 
unfaired roll-rate measurements, the differences between the motions in fig- 
ures l6(a) and l6(b) were attributed to the differences in roll-rate input data. 
The motions computed from faired roll- rate measurements are presented in fig- 
ure 16(c) . These motions compared favorably with analog simulation results 
which generated angular rates with marked simularity to the measured rates up 
to V55 seconds as illustrated in figure 18 . 

Final initial conditions .- By using faired roll-rate input data, trial- 
and-error integration of the flight measurements were made through maximum 
dynamic pressure until initial conditions at fifth- stage ignition were found 
that provided the best agreement with Bermuda optical tracking measurements and 
the wind-tunnel aerodynamic characteristics of the spacecraft. The final ini- 
tial conditions assumed for the analysis are presented in the following table: 


Table V 


Quantity 

Initial condition at - 

Ignition 

Burnout 

Flight time, sec 

377-85 

420.85 

V D , ft/sec 

17 761 

22 227 

ho, ft 

557 800 

308 082 

7 0 > de S 

-15.69 

- 20.25 

de g 

48.63 

51.12 

0 O , de 6 • 

- 54.70 

- 24.02 

Vo, deg 

81.70 

80.37 

0 O > de e 

65-24 

33-40 

Uq, ft/sec 

11 479 

19 953 

v Q , ft/sec 

-12 683 

-9854 

w Q , ft/sec 

5512 

5928 

Xgo' « 

311 693 

887 446 

Y go’ ft 

272 448 

891 156 

A 0 > deg 

-68.82 

-67.18 

Io, deg 

32.92 

31.26 
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It should be noted that some of the radar measurements at fifth- stage 
ignition were revised for better agreement with Bermuda optical tracking meas- 
urements. The revised initial conditions are listed and compared with radar 
measurements in the following table: 


Table VI 


Quantity 

Radar 

value 

Revised 

value 

Difference, 

percent 

V Q , ft /sec 

17 104 

17 761 

-5.7 

V ft 

565 582 

557 800 

1.0 

7 0 , tie g 

-15.69 

- 12.80 

6-5 

e 0 > de § 

48.85 

44.50 

8.9 


Although it was difficult to assess the errors in the radar measurements 
at the time of spacecraft separation and loss of radar tracking, the differ- 
ences in the initial velocity and altitude are believed to be within the accu- 
racy of the measurements. Since the radar measurements of initial latitude and 
longitude were assumed to be approximately correct, the largest changes required 
to obtain the best agreement with Bermuda optical tracking measurements were in 
the radar flight-path and heading angle as shown in figure 19. A comparison of 
the time histories of latitude and longitude is presented in figure 20. 

Although the revised initial flight -path angle in figure 19 appeared to be a 
reasonable extrapolation of earlier radar measurements, the difference in the 
initial heading angles was significant. However, a better match with Bermuda 
optical measurements was obtained with the revised heading angle as shown in 
figure 21. Although both radar and revised initial conditions gave good agree- 
ment in the elevation plane, a large change in initial heading angle was 
required to obtain closer agreement in the azimuth plane . 

The dynamic-pressure time histories obtained from both flight-data inte- 
gration and six-degree-of-f reedom simulation with wind-tunnel aerodynamics are 
compared in figure 22. Agreement between integration and simulation values for 
both radar and revised initial conditions was good during the buildup to maxi- 
mum dynamic pressure. Both values of longitudinal-force coefficient derived 
from integration results were within the scatter band of wind-tunnel test 
results. The differences after maximum dynamic pressure, however, indicated 
that the best agreement was obtained with revised fifth- stage initial condi- 
tions. In addition, the dynamic pressure computed from the measured accelera- 
tion and constant value of Cx = 0.575 was independent of initial conditions 

and agreed very well with both the integrated and simulation results. This 
comparison verified wind-tunnel predictions of constant longitudinal-force 
coefficient with Mach number up to angles of attack of about 20°. 

Trajectory initial condition study .- In order to determine whether there 
was any other combination of initial conditions which would result in the final 
dynamic-pressure time history of figure 22, a six-degree-of-freedom analog 
study of the effects of altitude, velocity, and flight-path angle on the space- 
craft acceleration time history was made starting at fifth-stage burnout. The 
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longitudinal-force coefficient of = 0-575? derived from the final dynamic- 

pressure time history, was assumed to be constant with angle of attack and Mach 
number for this investigation. The spacecraft motions, therefore, had little 
effect on the acceleration time history since the displaced coning motions at 
burnout precessed toward the velocity vector and eventually damped to angles of 
attack of less than 20° early in the flight . 

Families of reentry acceleration time histories were computed for a range 
of burnout altitudes, velocities, and flight-path angles. Two points on the 
acceleration time histories were assumed as necessary conditions for matching 
the flight measurements. These points consisted of the magnitude and time of 
occurrence of the maximum acceleration a^ and maximum rate of change of 

acceleration siL,max- For each change of initial conditions, then, the magni- 
tude of a^ max and a^ ma) , were determined at the same time that these quan- 
tities were measured during the flight. A summary grid of these results is 
presented in figure 25 , which shows the variation of a^ and ^ max with 

burnout flight-path angles for various combinations of altitude and velocity at 
fifth- stage burnout. Figure 25 indicated that only one combination of alti- 
tude, velocity, and flight-path angle would result in an acceleration time his- 
tory that would match the magnitude of both an ^ a L,max ^hat were 

measured during the flight. This combination agreed very well with those 
determined from trial-and-error integration of the flight measurements, and 
verified the uniqueness of the final initial conditions used to compute the 
spacecraft trajectory and motions. 

For other values of longitudinal-force coefficient within the scatter band 
of wind-tunnel test data, the velocity, altitude, and flight-path angle com- 
puted from the final initial conditions would be in error by approximately the 
amounts shown in figure 24. These errors represent a range of burnout flight 
conditions that would match the measured accelerations using wind-tunnel data 
and indicated that the final burnout velocity from integration of the flight 
measurements was within ±420 feet per second with an altitude and flight-path 
angle range of ±6200 feet and ±0.4°, respectively. It is believed that these 
increments also indicated the overall accuracy of the computed trajectory, 
since these errors would be incremental throughout the coasting trajectory. 

Data noise .- One of the problems encountered in integration of the flight 
measurements was the sudden shifts in the orientation of the computed motions, 
and divergence of the resultant angle- of- attack time history, as illustrated 
in figures 25(a) and 25(b). These results were caused by arbitrary fairing of 
telemeter measurements through noise periods which occurred at the flight 
times shown. New initial conditions were required, therefore, after each noise 
interval, and these conditions were determined by assuming that the motions 
were always around the velocity vector as indicated by earlier analog simula- 
tion studies. Table VII presents the initial conditions after each noise 
period required for integration of the flight measurements through maximum 
dynamic pressure. 
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Table VII 


Quantity 

Initial conditions 

Flight time, sec 

434.05 

447.45 

V Q , ft/sec 

22 334 

18 552 

ho, ft 

205 734 

105 436 

y 0 > deg 

-20.43 

- 20.67 

deg 

51. 14 

51.40 

9o, deg 

- 13 . 6 o 

- 20.70 

t 0 , deg 

67-35 

47.60 

0o> deg 

161.25 

231.56 

u Q , ft/sec 

21 890 

19 334 

v Q , ft /sec 

7473 

-29 

w Q , ft/sec 

-340 

526 


RESULTS AND DISCUSSION 


Reentry Trajectory 

Figure 2 6 presents the time histories of the spacecraft relative velocity 
and altitude computed from integration of the flight measurements using the 
final initial conditions of table V. These results are compared with measure- 
ments obtained from radar which was tracking the combined fourth and fifth 
stage before separation, and the burned-out fourth stage after separation. 
Comparison between radar and computed quantities shows that the spacecraft 
separated from the fourth stage and accelerated to a velocity of about 
22 230 feet per second at burnout of the fifth- stage motor at an altitude of 
about 308 000 feet. These results indicated that the fourth- stage motor fail- 
ure caused a loss in experimental velocity of about 3800 feet per second at the 
beginning of the heating data period near burnout at about 420 seconds flight 
time. The altitude at this time was higher than desired by approximately 
51 000 feet. The variations of the Reynolds number and dynamic pressure with 
Mach number obtained during the experimental heating data period are presented 
in figure 27 * 


Spacecraft Motions 

The motions of the spacecraft obtained from integration of the flight 
measurements during thrusting and coasting flight are shown by the Euler 
angle 9,\|/ cross plots in figures 28(a) to 28(d). The motions are plotted with 
respect to a fixed gravity-axis system since the rotation of the axis system 
traveling with the spacecraft was negligible during the time intervals shown. 
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Figure 28(a) shows that the displacement of the coning axis at separation 
of the spacecraft was about 25° below and 5° to the right of the initial veloc- 
ity vector. The coning half-angle at this time was approximately 26° and was 
damped during thrusting of the fifth- stage motor to about l4° at burnout. 
Damping of the motions did not occur symmetrically about a coning axis fixed in 
space, but was about an axis which precessed in the yaw direction and indicated 
the presence of asymmetric external moments during motor burning. 

Figures 28(b) to 28(d) show the residual coning motions at burnout which 
precessed toward and eventually damped around the velocity vector as the aero- 
dynamic restoring moments began to build up. It is interesting to note that 
the general shape of the motion pattern near burnout could have been estab- 
lished by integration of the flight measurements from the origin. For a given 
initial velocity vector, given attitude angle, and given attitude rate, the 
desired precession could then have been obtained by graphically translating and 
rotating the basic motion pattern as shown in figure 29. 

Figure 30 presents the time history of the spacecraft angles of attack 
associated with the motions shown in figure 28. These results indicated that 
the spacecraft began coning at separation with a maximum resultant angle of 
attack of about 49°. The angle of attack decreased to about 34° during the 
first 27 seconds of thrusting, and then began to build up to a maximum of about 
40° at burnout. This buildup appeared to start near the beginning of thrust 
tailoff and continued through burnout to about 426 seconds. 

A comparison of the envelopes of the resultant angle of attack after burn- 
out computed from the flight measurements and wind-tunnel aerodynamics is pre- 
sented in figure 31- During the time nose-cap calorimeter data were being 
recorded from burnout to about 438 seconds, the analog simulation results 
showed good agreement with flight data integration. The maximum angle of 
attack during this time varied from 40° to about 23° and indicated that the 
spacecraft was statically and dynamically stable with the nose cap on. After 
the nose cap burned off, the results from flight-data integration showed a sud- 
den change in amplitude of the oscillations which indicated that bumoff of 
the nose cap imparted an impulsive moment to the spacecraft in a direction to 
reduce the amplitude of the oscillations. Since the simulation did not account 
for this event, the angle-of -attack envelope continued at a higher level with 
amplitudes and trends similar to the integration results. Figures 30 and 31 
show that the spacecraft angles of attack after the nose cap burned off con- 
tinued to damp to a maximum of about 6° at maximum dynamic pressure. These 
results then indicated that the spacecraft was statically and dynamically 
stable throughout the experimental heating data period. 


Aerodynamic Results and Comparisons 

A check on the amplitude and phasing of the computed angle-of -attack time 
history, as well as the computed dynamic pressure, was obtained from comparison 
of wind-tunnel test results with force and moment coefficients derived from the 
accelerometer and rate-gyro flight measurements. Figure 32(a) presents the 
envelopes of the resultant-force -coefficient oscillations during flight in the 
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atmosphere with dynamic pressures greater than 100 lb/ sq ft. Resultant angle- 
of-attack envelopes predicted by wind-tunnel tests for the flight resultant- 
force coefficients are compared with the angles of attack computed from inte- 
gration of the flight data in figure 32(b). Fairly good agreement between the 
envelopes substantiated the trends and amplitudes of the final computed angle- 
of- attack time history. 

The phasing of the computed angle-of-attack oscillations was checked by 
comparison of the flight force and moment coefficients with the static wind- 
tunnel results. Figure 33 presents comparisons with the nose cap on and off 
through a Mach number range from about M = 20.6 to M = 9*7- The arrows on 
the continuous flight data show the manner in which the spacecraft computed 
forces and moments varied with time. For a known configuration immediately 
before and after the nose cap burned off at about M = 20.5; the flight force 
and moment variations showed fairly good agreement with static wind-tunnel 
results as shown in figures 33( a ) and 33(b) • Agreement at the minimum and max- 
imum oscillation peaks was also good. Although the flight data included 
dynamic effects, the similarity of the nonlinear characteristics of the data 
was especially noted. The resultant moment coefficients obtained from flight 
data in figure 33( a ) indicated the presence of a trim and hysteresis effect 
prior to bumoff of the nose cap. After the transient moments due to burnoff 
of the nose cap disappeared, agreement between the wind-tunnel and flight 
moments was exceptionally good as shown in figure 33(b). Figure 33(c) is pre- 
sented to illustrate the erratic variation of the computed forces and moments 
after the nose cap came off. Part of this behavior at the higher Mach numbers 
may be attributed to changes in the nose configuration as the result of abla- 
tion. Phasing errors or dynamic effects in the computed resultant angles of 
attack during this part of the flight may account for some of the discrepancies 
in the force comparisons. In general, however, the variations in the flight 
forces and moments were centered about the wind-tunnel static test results. 

Figure 3^ presents another comparison between the normal- and longitudinal- 
force coefficients from flight data and from wind-tunnel results. The flight 
normal-force coefficients and angles of attack in figure 3M a ) were obtained 
from oscillation envelopes for angle of sideslip 3 approximately equal to 
zero. Again, the wind-tunnel results represented a good average for the com- 
puted flight data. The longitudinal-force coefficients are compared in fig- 
ure 3Mb) and showed that the flight values were within the scatter band of the 
results from various test facilities. 


CONCLUDING REMARKS 


An analysis of the Scout ST-8 spacecraft reentry trajectory and angles of 
attack during the experimental data period has been made based on partial down- 
range optical tracking data, preflight wind-tunnel aerodynamic test results, 
and flight measurements of the spacecraft linear accelerations and angular 
velocities. The results indicated by this analysis are summarized as follows: 

1. The relative velocity attained by the spacecraft at the beginning of 
the primary heat-transfer data period was about 22 230 feet per second at an 
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altitude of approximately 308 000 feet. The initial experimental velocity was 
3800 feet per second lower than nominal at an altitude which was 51 000 feet 
higher as the result of fourth-stage motor malfunction. 

2. As the result of disturbances from fourth- stage motor failure, the 
spacecraft began coning at separation with a half-angle of about 26°. The 
coning axis at this time was displaced about 25° below and 5° to the right of 
the velocity vector. The maximum resultant angle of attack occurred near igni- 
tion and was approximately 49°. During the thrusting period, the maximum 
angles of attack damped to about 34° and then increased to about 4o° as the 
spacecraft entered the heating data period near burnout. 

3. The spacecraft maximum angles of attack during the primary heat- 
transfer data period varied from 4o° near burnout to about 23° at the time the 
nose-cap calorimeter burned off. After this time, the maximum angles of attack 
continued to damp through the ablation data period to about 6° at maximum 
dynamic pressure. These results indicated that the spacecraft was statically 
and dynamically stable throughout the heating data period. 

4. The force and moment coefficients corresponding to the computed tra- 
jectory and angles of attack were generally in good agreement with preflight 
wind-tunnel test results. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., January 28, 


1965 . 
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APPENDIX A 


EQUATIONS AND COMPUTATIONAL PROCEDURE FOR INTEGRATION 
OF RATE -GYRO AND ACCELEROMETER MEASUREMENTS 


As described in reference 2, body-axis equations of motion with a local 
gravity-axis system fixed at the spacecraft center of gravity (fig. 2) were 
programed for a spherical rotating earth on an IBM 7090 electronic data proc- 
essing machine. Rate-gyro and accelerometer flight measurements were continu- 
ously recorded on digital tape and used as input quantities in the equations of 
motion. 


The accelerometer measurements, corrected for errors caused by displace- 
ment of the instrumentation from the spacecraft center of gravity (ref. 2), 
provided the body-axis external forces (total force less gravity force) as 
illustrated by the following equations: 


So a L 


_ ?£ _ 


m 


m 


Fy Gy 

So 8 ! “ In Ur 

F7 G7 
g o a N = - “5T + it 


J 


(Al) 


where Gy, Gy, and G^ are the body-axis components of the spacecraft gravity 

force. The gravity forces were determined by introducing three unknown direc- 
tion cosines as follows: 


Gy = — W sin 9 
Gy = W cos 0 sin 0 
G^ = W cos 0 cos 0 




J 


(A2) 


With the use of equations (Al) and (A2) , the body-axis inertial velocities 
were obtained from integration of the flight measurements (p, q, r, ap, ap, 

and ajj) with the following basic equations: 
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u = u Q + 


(g 0 a L “ So sin 0 + vr - wq)dt 




V = V D + 
W = w D + 


/ 

/ 


(§o a T + So cos 0 sin 0 + wp - ur)dt 
(~So a N + So cos 0 cos 0 + uq - vp)dt 




J 


(A3) 


In order to determine the Euler angle direction cosines required for solu- 
tion of equations (A3) , the measured angular velocities p, q, and r were 
used in the following equations defining the spacecraft Euler angles: 


e = e 0 + r 

o 

* = to + / 

0 = 0 ° + J' 


(q cos 0 - r sin 0 + K@)dt 

( r cos 0 + q sin 
cos 9 

f sin 0 + K^dt 





J 


(A4) 


where Kg, and are the body-axis components of the angular velocity 

of the local gravity-axis system moving with the spacecraft with respect to an 
inertial reference. These quantities were a function of the spacecraft lati- 
tude and longitude as shown by the following equations: 


where 


Kg = L sin i|r - (co e + A) cos L cos i|r 

Ky = L cos ijr tan 0 + (<% + A) (sin L + cos L sin tan 0) \ (A5) 


H 


L + («e + *)cos L 

cos 9 cos 9 


J 


L = Lo + 
A = A 0 + 





(A6) 
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and 


L = - 


L g 

R q + h 


"\ 


A = 


X, 


g 


^R Q + h^cos L 


- ou 


(A7) 


Transformation of the body-axis inertial velocities to the local gravity- 
axis system was determined as follows: 


Xg 


Zq m-j_ nq 


u 

*g 

= 

l 2 m 2 n 2 


V 



J 3 

1 


w 


(a8) 


where Zq, mq, and nq are the Euler angle direction cosines between the 

spacecraft body-axis system and the local gravity-axis system as determined by 
the following equations: 


1 1 

1 2 

1 3 

m l 

HI2 

m3 

n l 

n 2 

n 3 


cos 0 cos 
cos 0 sin 
-sin 0 
sin 0 sin 
sin 0 sin 
cos 9 sin 
sin 0 cos 
sin 0 cos 
cos 0 cos 


* 

0 cos \jr - 
0 sin \(r + 
0 

0 cos \|r + 
0 sin \|/ - 
0 


cos 0 sin \[r 
cos 0 cos t 

sin 0 sin \|r 
sin 0 cos 


A 




J 


(A9) 
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and 


/*■ 


X g = X So + ./ X 6 dt 


Y g = Y So + J \ dt 


ZfT — Zcr _ + / Zrp dt 


J g ■ 


'g 


"\ 


(A10) 
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The spacecraft flight-path and heading angles were then computed from the 
following equations: 



X g - (Rq + h)tOe cos L 


"s 




J 


(All) 


The magnitude of the earth- relative velocity in equation (A12) was determined, 
the assumption being made that no winds were present through the altitude range 
of the experiment, from the following equations: 


V = |/u^ + v^ 


+ w 2 


(A12) 


where the earth-relative velocity components along the body axis were found 
from 


u = u - R 0 + h)o) e cos L 
v = v - m-j_(R 0 + h)<% cos L 
w = w - n-^(Ro + h)tD e cos L 


> 


(A13) 


The resultant angle of attack of the spacecraft was computed from the 
earth- relative velocity components in equation (A13) as follows: 




1 \/ v^ + w^ 

tan x 

u 


(Al4) 
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where the component angles of attack and yaw were defined by: 


a = tan - -'- ™ 
u 


p = tan -1 - ■ - — 

/ u 2 + w 2 




(A15) 


In order to compare the results from integration of the flight measure- 
ments with preflight wind-tunnel test results, the aerodynamic forces during 
coasting flight in the atmosphere were computed simultaneously from the fol- 
lowing equations: 


Wa 


c x - 


L,D 




qS 


(A16) 


c H = |/c B 2 + c Y 2 . ILf/ + ^ 




where the dynamic pressure q was determined by using the earth-relative 
velocity in equation (A12), and the free-stream density measured at the alti- 
tude computed from equations (A10) : 


q = | PV 2 (AIT) 

The thrust of the fifth- stage spherical rocket motor was determined from 
the following equation, negligible aerodynamic forces at the higher thrusting 
altitudes being assumed: 


T 


Wa L,T 


(Al8) 


The resultant moment coefficient during coasting flight in the atmosphere 
was found from the equation: 





(A19) 


It should be noted that the coefficients computed from the rate-gyro and 
accelerometer measurements in equations (Al6) and (A19) were the summation of 
the forces and moments including effects due to pitching and yawing angular 
velocities. 
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APPENDIX B 


EQUATIONS AND COMPUTATIONAL PROCEDURES FOR SIMULATION 
OF RATE-GYRO AND ACCELEROMETER FLIGHT MEASUREMENTS 


Six-degree-of -freedom equations of motion over a spherical rotating earth 
were programed on an arlalog computer using nonlinear aerodynamic character- 
istics of the spacecraft from preflight wind-tunnel tests of reference 3- 
Equations (A3) were used to compute the body-axis inertial velocities during 
coasting flight in the atmosphere by substitution of the spacecraft aerodynamic 
forces for the accelerometer measurements as follows: 


u = u Q + J ^ f x a “ ®o sin 9 + vr - wqjdt 

v = v Q + J ^Fy + g Q cos 6 sin 0 + wp - urjdt 

w = w Q + J ^Z A + ®o cos ® cos 0 + uq - vp^dt 


where the aerodynamic forces were computed from 


F X a = C X^ S 


~\ 


F y = CyqS = (-C R sin T + C N ™ qS 


F Z a = - C N^ S = (" C R cos T " % 


J 


> (Bl) 


The aerodynamic phase angle t was defined by the equations: 


(B2) 


i v _l w 

t = sin ^ — = cos - 1 - 


V v^ + w^ j v^ + w^ 


(B3) 


The angular velocities p, q, and r required for solution of equa- 
tions (Bl) , and comparison with rate-gyro measurements, were determined from 
the following equations written for a symmetrical spacecraft with no products 
of inertia and Iy = 1^: 
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where the total aerodynamic moments of the spacecraft were defined as 


% = 0 


M- 


■y a 


% ' 



(' 


sin t 



qSd 


qSd 


(B5) 




The force coefficients C x and Cp in equations (B2) and (B5 ) , and the 
spacecraft center of pressure x C p in equations (B5) were programed as a func- 
tion of Mach number and angle of attack from the wind-tunnel tests of refer- 
ence 3* The results obtained on both the model with the nose cap on, 

and the model ^N^Bi^ with nose cap off were used. The damping force and moment 

coefficients were estimated from Newtonian impact theory and were assumed to be 
equal in magnitude = Cy and C m ^ = Cn r ^ for the symmetrical spacecraft. 

The spacecraft center-of-gravity position in equation (B5) was obtained from 
preflight measurements with and without the nose cap installed. 

The spacecraft linear accelerations were computed for comparison with 
accelerometer flight measurements from equations (B2) as follows: 


% C x qS 
&L = ~ = ~W~ 


^ = CyqS 
a T W 


W 


% C N qS 

a N “ “ _ 

a W W 


*3 


(b6) 


The remaining equations defining the trajectory and motions of the space- 
craft were the same as those used for integration of the rate-gyro and accel- 
erometer measurements as presented in appendix A. 
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Figure 1.- Scout heat-transfer spacecraft. All dimensions are in inches unless otherwise noted. 










Figure 2.- Spacecraft coordinate systems and sign convention. 
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Figure 3* - Preflight nominal trajectory and flight events 
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Figure 5*- Continued. 



Figure 5.- Continued. 
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Figure 5.- Continued. 
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(a) t = 420-85 to 450.83 sec. 


(b) t = 430.85 to 435.83 sec. 



(c) t = 435.85 to 440.83 sec. 

Figure 7*- Cross plots of measured normal and transverse acceleration corrected for 

instrumentation displacement error. 
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(a) t = 440.85 to 1*52.83 sec. 
Figure 7*- Continued. 
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Figure 7.- Continued. 







Figure 7*- Continued. 






a^, g units 

(i) t = 472.85 to 477.85 sec. 


a T> g units 

(j) t = 477.85 to 482.85 sec. 



a T , g units 


(k) t = 482.85 to 488.49 sec. 
Figure 7-- Concluded. 
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Figure Variation of free-stream density with, altitude. 
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Figure 10.- Variation of measu: 



48 


Figure 11.- Comparison "between nominal and computed consumed weight during fifth- stage thrusting. 

(NOTS 100B) 17- inch spherical motor. 
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Figure 12.- Time histories of computed spacecraft weight, center of gravity, and moments 

of inertia during fifth- stage thrusting. 
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Figure 13-- Comparison 'between nominal and computed flight thrust time histories. 
(NOTS 100B) 17- inch spherical motor. 
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Figure 17.- Variation of computed resultant angle of attack 
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Figure 20.- Time histories of computed longitude and latitude. 
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Figure 22.- Comparison of computed dynamic -pres sure time histories. 
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(a) Euler angle displacement. 



(b) Resultant angle-of-attack displacement. 

Figure 25*- Shifts in Euler angles and resultant angle of attack during integration 

of flight data. 





Figure 26.- Time histories of computed velocity and altitude. 











o Jaqtunu sp-fouAey 

x — • 1 1 i 1 , 

'£> U\ J- CO CM r-i o 

qj bs/qp ‘aanssajd otuieuAq 


6k 


Figure 27.- Variation of computed Reynolds number and dynamic pressure with free- stream Mach number. 
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Figure 28.- Computed spacecraft motions. 
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Figure 31.- Comparison between computed resultant angle-of-attack envelopes 
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Figure 33 •- Continued. 
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(b) Comparison of longitudinal-force coefficient from M = 24.4 to 9 . 6 . 

Figure 34.- Comparison of normal- and longitudinal-force coefficients at oscillation peaks. 
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